Identification and characterisation of endogenous Avian Leukosis Virus subgroup E (ALVE) insertions in chicken whole genome sequencing data

Background: Endogenous retroviruses (ERVs) are the remnants of retroviral infections which can elicit prolonged genomic and immunological stress on their host organism. In chickens, endogenous Avian Leukosis Virus subgroup E (ALVE) expression has been associated with reductions in muscle growth rate and egg production, as well as providing the potential for novel recombinant viruses. However, ALVEs can remain in commercial stock due to their incomplete identification and association with desirable traits, such as ALVE21 and slow feathering. The availability of whole genome sequencing (WGS) data facilitates high-throughput identification and characterisation of these retroviral remnants. Results: We have developed obsERVer, a new bioinformatic ERV identification pipeline which can identify ALVEs in WGS data without further sequencing. With this pipeline, 20 ALVEs were identified across eight elite layer lines from Hy-Line International, including four novel integrations and characterisation of a fast feathered phenotypic revertant that still contained ALVE21. These bioinformatically detected sites were subsequently validated using new high-throughput KASP assays, which showed that obsERVer was highly precise and exhibited a 0% false discovery rate. A further fifty-seven diverse chicken WGS datasets were analysed for their ALVE content, identifying a total of 322 integration sites, over 80% of which were novel. Like exogenous ALV, ALVEs show site preference for proximity to protein-coding genes, but also exhibit signs of selection against deleterious integrations within genes. Conclusions: obsERVer is a highly precise and broadly applicable pipeline for identifying retroviral integrations in WGS data. ALVE identification in commercial layers has aided development of high-throughput diagnostic assays which will aid ALVE management, with the aim to eventually eradicate ALVEs from high performance lines. Analysis of non-commercial chicken datasets with obsERVer has revealed broad ALVE diversity and facilitates the study of the biological effects of these ERVs in wild and domesticated populations

A high-density SNP panel reveals extensive diversity, frequent recombination and multiple recombination hotspots within the chicken major histocompatibility complex B region between BG2 and CD1A1

Background: The major histocompatibility complex (MHC) is present within the genomes of all jawed vertebrates. MHC genes are especially important in regulating immune responses, but even after over 80 years of research on the MHC, much remains to be learned about how it influences adaptive and innate immune responses. In most species, the MHC is highly polymorphic and polygenic. Strong and highly reproducible associations are established for chicken MHC-B haplotypes in a number of infectious diseases. Here, we report (1) the development of a high-density SNP (single nucleotide polymorphism) panel for MHC-B typing that encompasses a 209,296 bp region in which 45 MHC-B genes are located, (2) how this panel was used to define chicken MHC-B haplotypes within a large number of lines/breeds and (3) the detection of recombinants which contributes to the observed diversity. Methods: A SNP panel was developed for the MHC-B region between the BG2 and CD1A1 genes. To construct this panel, each SNP was tested in end-point read assays on more than 7500 DNA samples obtained from inbred and commercially used egg-layer lines that carry known and novel MHC-B haplotypes. One hundred and one SNPs were selected for the panel. Additional breeds and experimentally-derived lines, including lines that carry MHC-B recombinant haplotypes, were then genotyped. Results: MHC-B haplotypes based on SNP genotyping were consistent with the MHC-B haplotypes that were assigned previously in experimental lines that carry B2, B5, B12, B13, B15, B19, B21, and B24 haplotypes. SNP genotyping resulted in the identification of 122 MHC-B haplotypes including a number of recombinant haplotypes, which indicate that crossing-over events at multiple locations within the region lead to the production of new MHC-B haplotypes. Furthermore, evidence of gene duplication and deletion was found. Conclusions: The chicken MHC-B region is highly polymorphic across the surveyed 209-kb region that contains 45 genes. Our results expand the number of identified haplotypes and provide insights into the contribution of recombination events to MHC-B diversity including the identification of recombination hotspots and an estimation of recombination frequency

CD99 and the Chicken Alloantigen D Blood System

The chicken D blood system is one of 13 alloantigen systems found on chicken red blood cells. Classical recombinant studies located the D blood system on chicken chromosome 1, but the candidate gene was unknown. Multiple resources were utilized to identify the chicken D system candidate gene, including genome sequence information from both research and elite egg production lines for which D system alloantigen alleles were reported, and DNA from both pedigree and non-pedigree samples with known D alleles. Genome-wide association analyses using a 600 K or a 54 K SNP chip plus DNA from independent samples identified a strong peak on chicken chromosome 1 at 125–131 Mb (GRCg6a). Cell surface expression and the presence of exonic non-synonymous SNP were used to identify the candidate gene. The chicken CD99 gene showed the co-segregation of SNP-defined haplotypes and serologically defined D blood system alleles. The CD99 protein mediates multiple cellular processes including leukocyte migration, T-cell adhesion, and transmembrane protein transport, affecting peripheral immune responses. The corresponding human gene is found syntenic to the pseudoautosomal region 1 of human X and Y chromosomes. Phylogenetic analyses show that CD99 has a paralog, XG, that arose by duplication in the last common ancestor of the amniotes

Genetic variation within the Mx gene of commercially selected chicken lines reveals multiple haplotypes, recombination and a protein under selection pressure.

The Mx protein is one of the best-characterized interferon-stimulated antiviral mediators. Mx homologs have been identified in most vertebrates examined; however, their location within the cell, their level of activity, and the viruses they inhibit vary widely. Recent studies have demonstrated multiple Mx alleles in chickens and some reports have suggested a specific variant (S631N) within exon 14 confers antiviral activity. In the current study, the complete genome of nine elite egg-layer type lines were sequenced and multiple variants of the Mx gene identified. Within the coding region and upstream putative promoter region 36 SNP variants were identified, producing a total of 12 unique haplotypes. Each elite line contained from one to four haplotypes, with many of these haplotypes being found in only one line. Observation of changes in haplotype frequency over generations, as well as recombination, suggested some unknown selection pressure on the Mx gene. Trait association analysis with either individual SNP or haplotypes showed a significant effect of Mx haplotype on several egg production related traits, and on mortality following Marek's disease virus challenge in some lines. Examination of the location of the various SNP within the protein suggests synonymous SNP tend to be found within structural or enzymatic regions of the protein, while non-synonymous SNP are located in less well defined regions. The putative resistance variant N631 was found in five of the 12 haplotypes with an overall frequency of 47% across the nine lines. Two Mx recombinants were identified within the elite populations, indicating that novel variation can arise and be maintained within intensively selected lines. Collectively, these results suggest the conflicting reports in the literature describing the impact of the different SNP on chicken Mx function may be due to the varying context of haplotypes present in the populations studied

Location of SNP within the promoter and untranslated region of chicken Mx.

<p>galGal4 Mx promoter region starting 400 nt upstream of the RNA transcription initiation site through the first 140 nt of the RNA transcript (excluding 4,647 nt from intron 1) are shown above. The previously described ISRE (−52 to −63) as well as other potential functional elements (ISRE2, −282 to −293; SP1-like element, −135 to −142; TATA box-like element, −19 to −12) are underlined. The “GAAA” motif found repeated in many IFN regulated gene promoter regions are shown in bold. SNP found in the 9 elite lines that differ from galGal4 are shown under the reference sequence. Additionally, the RNA transcription initiation and Mx 5′UTR, comprised of exon 1 and the first 92 nt of exon 2, is shown in italics.</p

Evidence of sites within the chicken Mx under selection.

a<p>chicken Mx codons with evidence of diversifying selection (dNS > dS) or purifying selection (dS > dNS).</p>b<p>p-value, *significance ≤0.1.</p>c<p>Bayes factor, * significance >50 significant.</p>d<p>Posterior probability, *significance ≥0.9.</p>e<p>not significant.</p><p>Evidence of sites within the chicken Mx under selection.</p

Location of SNP within the promoter and untranslated region of chicken Mx.

<p>galGal4 Mx promoter region starting 400 nt upstream of the RNA transcription initiation site through the first 140 nt of the RNA transcript (excluding 4,647 nt from intron 1) are shown above. The previously described ISRE (−52 to −63) as well as other potential functional elements (ISRE2, −282 to −293; SP1-like element, −135 to −142; TATA box-like element, −19 to −12) are underlined. The “GAAA” motif found repeated in many IFN regulated gene promoter regions are shown in bold. SNP found in the 9 elite lines that differ from galGal4 are shown under the reference sequence. Additionally, the RNA transcription initiation and Mx 5′UTR, comprised of exon 1 and the first 92 nt of exon 2, is shown in italics.</p